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On recent calculations of resonances for the Stark effect in hydrogen 
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We show that the resonances of the Stark effect in hydrogen reported in Phys. Rev. A 88, 022509 
(2013) are considerably less accurate than the number of digits appear to suggest. In particular, the 
imaginary part of the lowest resonance is several orders of magnitude greater than it should be. We 
compare the results of that paper with those provided by the Riccati-Pade method, perturbation 
theory and an asymptotic expansion for the resonance width. The inaccuracy of those results can be 
traced back to the lack of precision in the calculation of the matrix elements of the secular equation. 
We carry out a more accurate calculation with the same method and show that the agreement with 
earlier results of other authors is greatly improved. 


PACS numbers: 31.15.ac, 32.60.+i, 52.70.Ds 




In a recent paper Fernandez-Menchero and Summers [1| obtained the complex eigenvalues and eigenfunctions of 
the Hamiltonian operator for the hydrogen atom in a uniform electric field. They resorted to the Laguerre-mesh 
basis set proposed by Lin and Ho 2J for the treatment of the Yukawa potential in a uniform electric field and the 


complex-rotation method[3|. For brevity we will refer to this method as CRLM from now on. They compared their 

□ n n n 

results with those obtained by Lin and Ho[2J, Kolosov[4j, Rao and Li[5| and Ivanov|6| and overlooked the earlier 
impressive calculations of Benassi and GrecchiQ and the accurate results obtained by Fernandez jsj . Benassi and 
Grecchi resorted to a basis set of confluent hypergeometric functions that is suitable when the Sclrrodinger equation 
is written in squared parabolic coordinates. On the other hand, Fernandez applied the Riccati-Pade method (RPM) 
that does not require the use of complex coordinates. 

The results of Fernandez-Menchero and Summers [l] for the excited states deviate somewhat from the other ones m 
0] that appear to be in closer agreement among them. On the other hand, they were unable to compare their estimate 
of the lowest resonance for the field strength F = 0.005 a.u because there does not appear to be any independent 
calculation available. 

Fernandez-Menchero and Summers [l[ stressed that perturbation theory (PT) is unsound because the Stark effect 
does not fulfill the conditions for the application of such an approach. They are obviously unaware of the fact that 
the straightforward perturbation series is asymptotic to the real part of the resonances which in some cases can be 
obtained with reasonable accuracy by judicious truncation to a suitable number of terms. Besides, Adams [fj calculated 
the real and imaginary parts of the lowest resonance by means of analytic continuation and Pade summation and 
Jentschura[l0] resorted to Borel-Pade summation for obtaining several selected resonances. 

The purpose of this comment is to show that the resonances reported by those authors are considerably less accurate 
than the number of digits appear to indicate while the width for the lowest resonance is several orders of magnitude 
greater than it should be. The inaccuracy of their results, which accounts for the discrepancy with respect to those of 
other authors mentioned above, can be traced back to the lack of precision in the calculation of the matrix elements in 
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the secular equation. To verify this conjecture we carry out a more accurate calculation with the CRLM. In addition 
to it, we also resort to the RPM, PT and the asymptotic formula for the width of the lowest resonance derived by 
Benassi and GrecchiQ. 

The Schrodinger equation in atomic units is 


Hi/) = Eip 

H = -V 2 + --Fz , 

2 r 

where F is the intensity of the uniform electric field assumed to be directed along the z axis. 
This equation is separable in parabolic coordinates 

x = cos </>, y = yJTj] sin c/>, z = yy 

£ > 0, 77 > 0, 0 < <f> < 27r, 


(1) 


( 2 ) 


but the authors decided to treat the Schrodinger equation for the resulting Hamiltonian 


H = — 


Z + V 


d 
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d£ (5 <9£ ) + dr] dr] 


1 dr 
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as nonseparable. To this end they proposed the variational ansatz 
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^Nk(x) = (-!)■ y/Xk- 


( 3 ) 

( 4 ) 

( 5 ) 


X - Xk 

where Ln(x) is the Laguerre polynomial of degree TV and Xk its fcth zero. In order to obtain the resonances they 
resorted to the well known complex rotation method Q] that in this case is given by the particular complex-scaling 
transformation (£, rf) —> (e e l ®r]). The eigenvalues and expansion coefficients are given by the secular equation 

(H — ES)C = 0, (6) 

where the elements of the TV 2 x TV 2 matrices H and S are explicitly shown elsewhere^ and the elements of the column 
vector C are the coefficients Ckim■ Note that the integrals appearing in the matrix elements of both H and S should 
be calculated numerically and when we increase TV we have to calculate all those integrals again. 

On the other hand, in order to apply the RPM we write 


</’(£, V, 4>) = i€v) V rn = 0, ±1, ±2,..., 


( 7 ) 


and obtain two equations of the form 


cP_ 

dx 2 


1 — m 2 
Ax 2 


E F A,j 1 . . 

2 ~ a ^ x + ~ = °> 


( 8 ) 


where a = ±1 and A + = A and A_ = 1 — A are separation constants. When er = 1, x = £ and $(£) = u(£); when 
a = —1, x = 7] and $( 77 ) = v{rf). 
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The regularized logarithmic derivative 


ffa) 


can be expanded in a Taylor series 


s $'(a;) \m\ +1 

x d)(a:) ’ S 2 


(9) 


f(x) = fa 

j =0 


( 10 ) 


where the coefficients fj are polynomial functions of E and A. The details of the method are outlined elsewhere 
here it suffices to say that we construct the Hankel determinant 


where D = 2,3,. 
equations 



fd -\-1 

fd-\- 2 

fD+d 

H d D (E,A,F) = 

fd+2 

fd- 1-3 

• f D+d +1 


fD+d 

f D+ d + l 

• f2D+d-l 


and d = 0,1,... is kept fixed. The approximate eigenvalues are given by the set of nonlinear 


H d (E,A,F) = H d ) (E,A-l,F) = 0. (12) 

The main advantage of the RPM is the enormous rate of convergence of the approximate eigenvalues E^ D,d ^ which 
enables us to obtain very accurate results with determinants of relatively small dimension. However, the great number 
of roots in the neighborhood of each eigenvalue makes it difficult to find the optimal sequence that converges to it. 
We can mention two other disadvantages of the RPM: first, it only applies to separable problems; second, it does not 
give the eigenfunction but its logarithmic derivative. Consequently, it is not practical for the calculation of physical 
properties that are defined in terms of the eigenfunctions. In spite of these disadvantages the RPM is an extremely 
useful tool to produce accurate results for benchmark purposes in those cases where it can be applied. The Stark 
effect in hydrogen is one such problem as shown in what follows. 

In order to test the rate of convergence we calculate log | — of- D,d ^) /a^- D+1 ’ d ^ | for a = R eE or a = IvaE. 

Here we restrict ourselves to d = 0. The results for the lowest resonance when F = 0.005 shown in Fig. |T] strongly 
suggest that present RPM calculation may be remarkably accurate. 

The converged RPM resonance, truncated to a reasonable number of digits, is shown in Table Q] together with the 
CRLM one)l]. The discrepancy between the real parts (1 in 10 6 ) may not appear to be so serious at first sight, but the 
imaginary parts differ in many orders of magnitude predicting considerably different lifetimes for the same metastable 
state. Our estimated value of Im E is in perfect agreement with the analytic asymptotic formula derived by Benassi 
and Grecchijl]: 

ImE -2 F- 1 e~ 2/(3F) (1 - 8.916F + 25.57 F 2 + 0{F 3 )) , (13) 

also shown in Table [I] From a practical point of view this state may be considered stable because there is no experi¬ 
mental way of determining such a long lifetime. However, it is worth stressing the fact that the RPM enables us to 
calculate it. 
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In principle we expect that a properly truncated perturbation series will exhibit an accuracy of the order of 
|Im£'|. On summing the first 130 terms of the perturbation series calculated by means of the hypervirial perturbative 
method 11] we obtained the following result: 


E pt = —0.5000562847937929693317739476914328819632509273188913726, 


ReU 


RPM 


= -0.50005628479379296933177394769143288196325092731889137262135731287257, 


(14) 


that agrees with the RPM one to the last digit. Obviously, this result is more accurate than the CRLM one and most 
probably than any other calculation based on diagonalization and complex rotation. 

In our opinion it is quite difficult to obtain so small imaginary parts of resonances by means of standard methods 
based on diagonalization and complex rotation which may probably be the reason why it was not included in the 
earlier calculations chosen for comparison^, Q-fg| . This fact clearly shows that the RPM is a most valuable tool for 
testing other calculations on separable problems. 

The discrepancy between the results of Fernandez-Menchero and Summers and those chosen for comparison is 
probably due to insufficient accuracy in the calculation of the CRLM integrals. In order to test this conjecture we 
calculated those integrals with greater accuracy (15 digits) by means of a computer algebra system that enables, in 
principle, unlimited precision and then solved the CRLM secular equation in the usual way for N = 30. Present 
CRLM lowest resonance, also shown in Table (T) is in perfect agreement with the RPM result. Tabic |TT] shows that 


present CRLM resonances are in better agreement with those of the other authors 


(for brevity we only show 


the results of Lin and HoQ). In passing we mention that the width of the state | 2 1 0 ) calculated by the other 
authors^, 4- 61 was incorrectly transcribed by Fernandez-Menchero and Summers. 

Typically the minimum value of |Im.E| that one can calculate by a diagonalization method is of the order of the 
accuracy of the eigenvalue. The accuracy of the results of Fernandez-Menchero and Summers is roughly of the order 
of 10 -6 which is the reason why they estimated such width for the lowest resonance. If we take into account that 
other physical effects (magnetic-field perturbation, fine structure, etc.) already mentioned by the authors are of the 
order of 10~ 4 we conclude that the lack of precision just mentioned may probably not be serious for most practical 
applications. However, we have clearly shown that the CRLM results can be considerably improved by increasing the 
accuracy of the integrals that provide the elements of the matrices in the secular equation. 

We add that PT and the asymptotic formula ^ are remarkably useful for the calculation of extremely sharp 
resonances. Unfortunately, while we have the perturbation series for all the resonances [llj the asymptotic formula is 


only available for the lowest oneQ. It should be mentioned that, even for the small field intensity considered here, 


PT is not so accurate for higher resonances unless one resorts to a suitable summation method|9. llOl]. On the other 
hand, the RPM yields eigenvalues of similar accuracy to the one shown here. 
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TABLE I: Lowest resonance for F = 0.005 


ReE 

CRLM |4J -0.5000553416 

Present CRLM -0.500056284793 
RPM -0.5000562847938 

Asymptotic 


r 

0.8944475605 x 10~ 7 
< 1 x 10~ 13 
9.49802741674 x 10 -56 
9.4983 x 10" 56 


TABLE II: Resonances for F = 0.005 



ReE 

r 

|2-10) 



CRLM [1] 

-0.1426203564 

1.057292433 x 10” 4 

Present CRLM 

-0.1426186075727079 

1.05944463673 x 10“ 4 

Ref. [2] 

-0.1426186076 

1.059444711 x 10” 4 

1210) 



CRLM [1] 

-0.1120633027 

4.930560122 x 10~ 6 

Present CRLM 

-0.1120619240019936 

5.72936843930 x 10“ 6 

Ref. [2] 

-0.1120619240 

5.72939466 x 10' 6 

1 2 0 1 ) 



CRLM [1] 

-0.1271464039 

2.671348551 x 10“ 5 

Present CLRM 

-0.127146612703972 

2.6152854466430 x 10' 5 

Ref. [4] 

-0.127146612 

2.61528545 x 10 -5 



FIG. 1: Convergence of the RPM for the lowest resonance when F = 0.005 







